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Abstract 

Helium atoms or hydrogen molecules are believed to be strongly bound within the interstitial 
channels (between three carbon nanotubes) within a bundle of many nanotubes. The effects on 
adsorption of a nonuniform distribution of tubes are evaluated. The energy of a single particle 
state is the sum of a discrete transverse energy Et (that depends on the radii of neighboring tubes) 
and a quasicontinuous energy E z of relatively free motion parallel to the axis of the tubes. At low 
temperature, the particles occupy the lowest energy states, the focus of this study. The transverse 
energy attains a global minimum value (Et = E m i n ) for radii near 

Rmin = 9.95 A for H 2 and 8.48 

A for 4 He. The density of states N(E) near the lowest energy is found to vary linearly above this 
threshold value, i.e. N(E) is proportional to (E — E m i n ). As a result, there occurs a Bose-Einstein 
condensation of the molecules into the channel with the lowest transverse energy. The transition is 
characterized approximately as that of a four dimensional gas, neglecting the interactions between 
the adsorbed particles. The phenomenon is observable, in principle, from a singular heat capacity. 
The existence of this transition depends on the sample having a relatively broad distribution of 
radii values that include some near R m in- 



I. INTRODUCTION 



One of the most persistently interesting topics in condensed matter physics is Bose- 
Einstein condensation (BEC). First postulated |l| on the basis of an idealized model (non- 
interacting particles), BEC was subsequently argued by London to provide an explanation 
of superfluidity in 4 He, which is a system consisting of relatively strongly interacting atoms 
[2j]. The London hypothesis has been confirmed in both detailed calculations and careful 

n n 

neutron scattering experiments [3J, |4|. More recently, BEC has been observed to occur in 
ultracold gases consisting of alkali or H atoms and in systems of excitons at low temperature 

A logical candidate for the occurrence of BEC is para-H 2 , which is the boson equilibrium 
form of hydrogen at low T. The apparent nonexistence of BEC for bulk H2 is attributed to 
its crystallization below 14 K, which preempts BEC. Some three decades ago, Ginzburg and 
Sobyanin proposed a search for BEC and superfluidity of Ho in confining geometries (for 
which crystallization occurs at a lower temperature, if at all) £ L This su gge stion stimulated 
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a number of searches involving both theory and experiments |9|, 
There has been at least one assertion of successful observation of BEC, involving the lack of 
dampingof rotational motion of an impurity surrounded by hydrogen fluid within a small 
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In this paper, we describe a different geometry that is predicted to allow BEC to occur 

is an ensemble of carbon nanotube bundles. A crucial role 
is played in this phenomenon by the presence of nonuniformity. In particular, we assume 
that there exists a fairly broad distribution of nanotube radii present in the experimental 
sample. We focus our calculations on the spectrum of particles (either H2 or 4 He) confined 
within interstitial channels (ICs) formed by groups of three nanotubes. It has been argued 
that both of these species, if allowed access to such channels, are more strongly bound 
there than on the surface of graphite, the most strongly physisorbing planar surface (lflj . 
There is ambiguous experimental evidence concerning the relative binding energy of these 



species in the IC and in the grooves on the outside surface of the bundle |2C 
We interpret the very large isotope effect seen in isosteric heat data of Wilson et al [22j 
to indicate particle localization within the ICs for their sample. The specific behavior is 
evidently sample- dependent, just as our calculations are very sensitive to the assumptions 
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that we make. 

Much recent research has been devoted to hydrogen within nanotubes bundles, partly due 
;o the exciting prospects of technologies related to hydrogen storage and isotope separation 
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291 ] . The vast majority of the theoretical research to date (including that 



of our group) has assumed a monodisperse distribution of tubes. The present study was 
stimulated by recent results of Shi and Johnson jscj, who showed that the distribution of 
tube sizes present in actual samples leads to predictions that are quite different from those 
based on the monodisperse models, in better agreement with experimental data than the 
naive models' predictions. 

The outline of this paper is the following. Section II describes our calculations of the 
energy spectra of 4 He and H 2 gases within ICs. Section III employs these results to compute 
a set of predictions related to BEC that are testable in principle. Section IV presents a 
further discussion of the model and draws conclusions. 



II. ENERGY SPECTRA AND DENSITY OF STATES OF ADSORBED GASES 

We focus here on those states having the lowest energy since they are the most important 
for the BEC phenomenon; higher energy behavior can be explored with the same approach 
but minor changes in the analysis are required. The energy of a particle confined within an 
IC is evaluated with the Schrodinger equation. Assuming that the particular IC is perfectly 
straight, this equation is separable; then the total energy is given by the relation: 

E n (p) = E*' ans + E lon °(p) (1) 



E lon9 (p) = 2- 



(2) 



Here, the first term is the n-th eigenvalue of the Schrodinger equation describing motion 
transverse to the axis of the IC and the second term is the quasicontinuous kinetic energy of 
motion parallel to the axis. We assume for simplicity that the longitudinal kinetic energy is 
determined by the free particle mass; if the potential were corrugated the band mass should 



be used instead of m to characterize the low energy states 
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Such a substitution 



would lead to a straightforward change in the numerical and analytical results below. 



In the present circumstance, the transverse eigenvalues are separated by energies (of order 
100 K) that are much larger than relevant temperature scales |33j]. Hence, we need consider 
only the lowest eigenvalue E^ ans of transverse motion. We denote this quantity E t (R,), where 
the vector R = i? 2 , -R3) has components equal to the radii of the tubes surrounding the 
IC. We refer to the domain of possible ICs as R space; a given IC is represented by one 
point in this space while a sample with many ICs is described by a cloud of points in this 
space. The single particle spectrum is now specified by the relation 

E(R,p) = E t (R) + ^ (3) 

Here, E t (R) is computed following the procedure of Stan et.al.[34]. To get a sufficiently 
accurate dependence of E t on R, one must include anisotropic and anharmonic contributions 
to the potential in the IC. We note that alternative parametrizations of the potential have 
been proposed. Table 1 presents one such alternative result for E t . The qualitative behavior 
presented in the remainder of this paper is the same for the case of that alternative potential. 

Of considerable interest is the density of states h(E, R) for a single IC. For a specific 
set of adjacent tubes (R), this function is given by a sum over all momenta, which may be 
replaced by an integral: 

h(E, R) = X*E ~ E(R, p)] = ^) l ' \ E _ E ] ml/ , 8[I5 - MR)] (4) 

Here L is the length of the channel and the usual procedure for quantizing the ID motion 
has been followed. The inverse square root dependence on the energy above threshold is 
characteristic of ID motion. The last factor involves the Heaviside step function which 
is unity here for E > E t (R) and zero otherwise. 

Now let us consider an ensemble of nanotubes. The lowest energy states (those domi- 
nating the low temperature behavior) are concentrated near the minimum of the function 
E t (R). As might be expected, the global minimum of E t (Ii) is found for the symmetric case, 
{Ri} = Rmin- We define the vector specifying this IC as R m i„ = R m i n (l, 1,1) and E min as 
the corresponding energy. The numerical results for these quantities appear in Table 1. 

The behavior of the function E t (R) near this minimum is remarkable. If one considers 
only symmetric ICs (bounded by tubes of identical radii), the energy varies rapidly as a 
function of the difference |R — R m j n |. However, the spectrum of the very lowest energy 
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states is dominated by asymmetric ICs for which one of the radii equals R m in and the two 
others differ from R m in by equal, but opposite amounts, e.g. R = R m i n (l, 1 + x, 1 — x), 
where i<1. The energy varies extremely slowly with x in such a case. This is depicted for 
H 2 in Fig. 1, which shows the energy variation in a plane within R space that contains both 
this variable x line and the diagonal (1, 1, 1) line. (Analogous behavior occurs for the case of 
4 He). One observes in the figure a very narrow valley of low- lying states of variable x, along 
the (0, 1, —1) direction, where the origin is shifted to R m i n . Similar small gradient behavior 
of E t (FL) is found along equivalent permutations of this direction. Note that the higher 
energy contours are extended, nearly straight lines perpendicular to the diagonal. These 
latter contours imply an energy dependence of the density of states that is qualitatively 
different from that based on contours of lower energy (near the valleys). 

We next evaluate the transverse density of states for a given experimental sample of 
tubes. We denote by ^(R) the density distribution of ICs, defined so that i/(R)<£R is the 
number of ICs present in the sample within a volume dH in R space. In this paper, we 
focus on low energy behavior corresponding to such close proximity to R m j„ that we may 
replace z/(R) by z/(R m j n ) ; future studies will address more general situations for which 
that approximation is not adequate. In the following discussion, we assume that z^(R) is a 
sufficiently broad distribution that ICs near this minimum exist in sufficient numbers to treat 
adsorption statistically. Moreover, we assume that the set of tube sizes forms a continuum. 
The former assumption is well justified in experimental samples produced to date involving 
arge numbers of tubes, which typically have a dispersion in radius values of order 20% 



radius 1 nm, there are many values of the chiral indices that yield quite similar radii; the 
mean spacing between successive values of the radius is quite small, 0.002 nm. In addition, 
the ICs experience small perturbations due to their environment; in particular, the ICs at 
the center of a bundle are expected to be compressed relative to those near the perimeter. A 
model calculation (unpublished) yields an expression for the difference 5 between the tube- 
tube separation in an infinite rope and that {Ri) of an isolated nanotube pair; apart from a 
constant of order one, the result is 5 = [CA 2 ]/[/cRf], where C is the van der Waals-London 
intercarbon interaction coefficient [C ~ 20 eV-A 6 ), A is the ID density of carbon atoms in 
a nanotube (A ~15A _1 ), Ri ~17 A is the spacing between tube centers and fc=10~ 3 eV/ A 3 
is the force constant per unit length associated with neighboring tubes' interactions [37]. 




The latter (continuum) assumption is justified by the following argument. Near 
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The result, S ~0.15A, represents a nearly 5% compression of the lattice constant of ropium 
relative to the separation of an isolated pair. This shift in se par ation is consistent with 
estimates of the analogous shift in breathing mode frequency [38(. On the basis of these 
arguments, the assumption of a continuum of possible radii appears appropriate. 

The transverse density of states g(E) for a given sample of tubes is expressed by the 
relation 

g(E) = J dRv(R)5[E-E t (R)\ (5) 

The total density of states N(E) is found by summing over all states of the ICs present in 
the sample: 

N(E) = [ dRj2$[E - E(R,p)}v(R) ( 6 ) 
J v 



N(E) = 4(y) 1/2 l E dE t g(E t ) [E - E t }~ 1/2 (7) 
Tin / j E m i n 

Note that in the special case of a perfectly uniform distribution of Nic identical ICs, f(R) = 
N IC S(R - Rq) , yielding g(E) = N IC S[E - E t (R )]; then N(E) is precisely N IC times the 
single tubes density of states, h(E,R ), given in Eq. 4 above, as expected. 

Figs. 2 and 3 present the transverse densities of states g(E t ) for He and H2, respectively, 
computed numerically from a distribution of 600 3 tubes with radius values spread uniformly 
over the interval 8.5 A to 11.5 A. In both cases, the behavior of g(E) at low energy is 
characterized by a power law dependence on energy above threshold : gi ow {E) ~ (E — 
Emin) 1 ^ 2 ■ At higher energy, there is a decrease of g(E) with energy, fit well by the expression 
ghigh(E) ~ (E — Ernin)^ 1 /' 2 . The initial behavior is explained by the following argument. 
At low energy, the density of states is obtained by counting the number of ICs whose 
energy E t (R) lies in an interval between E and E + dE. Assume for simplicity that the 
isoenergy surfaces are spherical surfaces centered about the point R = R m i n (i.e. E — E min = 
(k/2)\R — R m j„| 2 ). Then g(E) is obtained from the number of points in the shell between 
spherical surfaces (corresponding to E and E + dE) centered at R m in- The result of this 
simple model is then 
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giow(E) = z/(R mm ) 4tt (p) 1 / 2 (E - E min ) 1 ' 2 (8) 

This dependence on the square root of the energy turns out to be valid even in the extreme 
anisotropic case of interest here j^. The result then is that k 3 is replaced by k 3 nis = k x kf, 
where x denotes the (1,1,1) direction in R space and t denotes the two directions transverse 
to that. 

Note that this square root energy dependence of gi ow (E) coincides with that of the density 
of states of a 3D gas. The higher energy behavior (1/y/E) of ghi g h{E) is quite different. The 
relevant high energy regime in Fig. 1 is that where the energy surfaces in R space are 
perpendicular to the (1, 1, 1) diagonal. Since this variation is essentially a ID dependence, 
the resulting high energy form of ghi g h{E) is that of a ID system, as in Eq. 4 above; this gives 
rise to a behavior ghi g h(E) ~ (E — E min )~ x l 2 (as is observed in Figs. 2 and 3). We remark 
that this energy-dependent variation in effective dimensionality has analogs in electronic 
band structures, where wave vector replaces the R variable as the source of the unusual 
dependence. |4I| 

Figs. 2 and 3, lower panels, depict the total densities of states N(E) derived from Eq. 
[7| In view of the power law forms (e.g. Eq. |HJ) of g(E), we expect power laws for N(E) in 
appropriate regimes of E. For example, the very lowest energy behavior is determined by 
the integral 

f E dE t {E t -E min ) l l 2 (E-E t )- 1 ' 2 = {E-E min ) f dy [y/(l - y)} 1 ' 2 = - (E - E min ) (9) 



N(E) = u(R min ) — (pH V2 ( E - E mvn) (10) 

"'anis 

Note that the prefactor of the energy difference (which we call a) is proportional to the 
volume of the nanotubes, so that N(E)/a is an intensive variable. The linear dependence 
on energy above threshold (E — E m i n ) is the behavior characteristic of a 4D gas. This 
initially unexpected result is a logical consequence of convoluting a transverse spectrum 
characteristic of a 3D system (eq. ISJ) with the ID longitudinal degree of freedom. Similarly, 
the behavior at higher energy is obtained by convoluting an inverse square root dependence 
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of ghigh(E) with the ID inverse square root dependence. Then the high energy result is 
that N(E) is independent of energy, behavior characteristic of a 2D density of states. These 
analytic results for the spectra are consistent with the numerical results seen in Figs. 2 and 
3; that is, the low energy density of states is proportional to (E — E min ) while the high 
energy behavior is independent of energy. By the word "high" , we mean energies reasonably 
close to E min (say E — E min ^20mK) but not the absolutely lowest energies. 

To summarize this section, we note that unexpected behavior of N(E) emerges from a 
simple model. The low energy spectrum of the nonuniform system is qualitatively different 
from the spectrum obtained when heterogeneity is ignored. The latter corresponds to ID 
physics, i.e. N(E) proportional to (E — E min )~ l l 2 . With heterogeneity, instead, we find that 
N(E) exhibits a 4D form at very low energy and a 2D form at higher energy. 



III. PREDICTION OF BEC 



In this section we explore the consequences of the anomalous densities of states for thermal 
properties of these systems. In so doing, we assume that the adsorbed gas can equilibrate. 
This entails the rearrangement of the particles as the temperature T is lowered by moving 
from ICs that have high transverse energies to those with lower energies. Thus, we assume a 
combination of sufficient particle mobility and sufficiently patient experimentalists that the 
assumption is valid. In the extreme opposite limit of negligible transport of particles the 
thermal behavior is that of isolated and independent ID systems. The difference between 
these behaviors is dramatic. 

At very low density, one can ignore the effects of both interparticle interactions and 
quantum statistics. Then, the specific heat C (T) = dU / dT at fixed number (N p ) of particles 
can be computed from the classical Boltzmann expression for the mean energy U: 

MT) JdEEN(E) ex P [-pE] 

N p > class%cal JdEN(E)}exp[-pE] 1 j 

Here = /c^T. One can insert possible power law forms N(E) ~ i^W 2 - 1 ) corresponding 

to a gas in d dimensions and derive the familiar relation 

/classical 7T \ I 



N p k B 2 
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Fig. 4 presents the resulting classical behavior for the systems of He and H2 in nanotubes. 
One observes that the low T regime does correspond to the expected 4D gas behavior, 
C/(N p ks) = 2 while the high T regime exhibits 2D gas behavior, C/{N p k,B) = 1- The initial 
rise above 2 at very low T is explained by a small upward deviation from linearity found 
in N(E) at low E. That is, if N(E) at low E is proportional to E(l + E/E\), where E\ 
is some constant, then [C / (N p kB)]ciassicai = 2[1 + 2/{J3E\)\. The crossover between the 4D 
and 2D regimes is gradual, beginning at a few mK and complete by 20 mK; the latter value 
corresponds to the energy above which N(E) is virtually identical to the 2D form. 

At higher density, one must consider the effects of interactions and quantum statistics. 
In the present treatment we ignore interparticle interactions, except for some comments in 
the next section. While this is a common assumption, usually adopted for simplicity, there 
exists some justification for it in the present instance because of electrodynamic screening 
of the van der Waals interaction by the surrounding medium; there are additional elastic 
screening effects that have yet to be studied 4.1]. We do not know the resulting effective 
interaction well enough to characterize it. In contrast, the effects of quantum statistics are 
easily evaluated by the standard procedures. In particular, the assumed constraint of fixed 
particle number yields an implicit relation for the chemical potential fi, which depends on 
T, from which one may derive the energy: 

N p = J dEn(E, T) (13) 
U(T) = J dE En(E, T) (14) 

n(E, T) = N(E) (exp[i3(E - //)] - l)" 1 (15) 

(As usual, the classical equations ^2 and are obtained in the limit when the fugacity 
Z = exp (/?//) <C 1). From U(T), one may compute C(T). In evaluating these formulae, 
one encounters a phenomenon that is well known to occur in 3D bose gases: BEC. This 
transition is derived here in the usual way, by evaluating the implication for the variable 
ji (at fixed N p ) as T is lowered: upon decreasing T, /1 increases up to the point where it 
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reaches its highest possible value /x = E min . Thus, the BEC condition for the transition 
temperature T c = (ksPc)^ 1 is 

N p = J dEN(E) (exp[(3 c (E - E mm )] - l)' 1 (16) 

Since \x cannot exceed E mini as T falls below T c (given implicitly by this relation), there 
ensues a division into two groups of particles- the Nexcited particles accounted for by the 
integral over excited states (E > E min ) and the remaining particles having E = E min that 
comprise the bose condensate, N cond : 

N cond = N p - N excited (T<T C ) (17) 

N excited = J dEN(E)(ex P [P(E - E mm )} - l)' 1 (18) 

The results of these calculations appear in Figs. 5 to 8. Fig. 5 depicts the dependence of Z 
on the relative temperature T/T c for two different densities. Fig. 6 displays the dependence 
on T of the condensate fraction f(T) = N con d/N p . Fig. 7 shows the function C(T), which 
exhibits singular behavior at the transition. Finally, Fig. 8 shows the dependence of T c on 
N p . Note that the transition temperature falls within the mK range studied in many low 
temperature laboratories, so experimental observation of this transition is feasible. 

In the preceding three figures, comparison is made with results obtained for a purely 4D 
bose gas defined to have the same initial slope of the density of states as the real system. 
That ideal system satisfies N(E) = aE, for all E > 0, with the analysis presented in the 
Appendix. (For the real system, the a coefficient is the prefactor of the energy term in Eq. 
HD|) . The key results for the 4D system are 

Ud{T) = 1 - {^-f (19) 

J-c,AD 

Cad{T) = ^(^) 2 , T < T cAD (20) 

T cAD = (^)V3_L (21) 
a ttkb 
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Here £(3) = 1.202 is a Riemann zeta function. The T 2 dependences of f(T) and C(T) 
are exhibited by the "real" system at low T, as seen in the figures. At the very lowest 
temperatures, only the lowest energy states play a role; since N(E) is the same for the ideal 
system in this region, the behaviors coincide. The transition temperature is higher for the 
real system than for the ideal one because at high E, N(E) is lower in the real case than in 
the ideal case, implying that (for the same value of N p ) fj, reaches its transition value E min 
at a higher value of T c . 

Since the ground state energy of the noninteracting classical system coincides with that 
of the noninteracting quantum system, as does the energy at very high T, it follows that the 
integral of the two heat capacities is the same. At any particular T, of course, the values of 
C are very different. In particular, the quantum system manifests the reduced C at low T 
expected in this bose gas case. The quantum system exhibits a much larger value of C in 
the vicinity of T c , which is dominated by a large rate of excitation out of the condensate; 
see Fig. 7. 

One interesting question is this: how are the particles distributed in energy as a function 
of T. To address this, we present in Fig. 9 the occupied energy distribution function n(E, T) 
at various temperatures. One observes the increasing relative occupation of low energy states 
with decreasing T; for T > T c , their sum is just N p , the total number present. At T c there 
occurs a discontinuous jump in the occupation of the lowest energy states accompanying the 
BEC transition. Below T c , the number N excite d of uncondensed particles decreases, vanishing 
at T = 0, when all of the particles settle into the condensate. 

Our prediction of BEC in this anisotropic system raises a question: what is the relation- 
ship between the transition condition and that of more familiar homogeneus problems, e.g. 
the 3D Bose gas and the Kosterlitz-Thouless transition of a superfluid film? In the latter 
examples, the relevant criterion is if the form p\ d ~ 1, where p is the number density, d the 
dimensionality and A = (2irf3h 2 /m) 1 / 2 is the de Broglie thermal wavelength. We present now 
a heuristic argument that is consistent with this homogeneus relation and yields a good es- 
timate of the transition temperature for the present BEC problem. The argument estimates 
T c by analyzing the condition for quantum degeneracy; the number of particles N p becomes 
comparable to the number of states N states (T c ) lying within fcfiT c of the ground state. The 
number of states within that interval is essentially the product of the number of states of 
longitudinal motion within that range and the number of states of transverse motion within 
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that range: 



N states (T) ~ iV||(T) JVj.(T) 

The factor jVn is the ratio of the thermal momentum in the z direction to the spacing 
between discrete states. The factor N± is derived from a thermal spread of states in R space 
[SR ~ {{3k an i s ) _1 / 2 ] and the corresponding number of states in the experimental sample, 
N± ~ uSR 3 . Setting N p = N st ates{T c ) yields our estimate of the BEC transition condition: 

Lv m x l 2 , . 

~ fofi2~jj]2~ (22) 

' '"ams 

Apart from numerical factors, this result coincides with the combination of Eqs. 10 and 21. 
Note that this last equation can be rewritten in the form p*X A ~ 1 with p* = N p /(L I 3 ) 
where I = v l l 3 Ti/{m k anis ) l l 2 is a characteristic length defined by the energetic parameters 
and the distribution of radii. 



IV. SUMMARY AND DISCUSSION 



In this paper we describe a conceptually simple system that turns out to exhibit remark- 
able behavior. The focus of the study is the lowest-lying states of a system of hydrogen 
or helium moving within a collection of nanotubes having many ICs. Our analysis yields a 
density of states N(E) that fulfills the requirement of BEC (i.e. the maximum integrated 
occupation number is finite). Equally interesting, perhaps, is the finding that the calcu- 
lated properties at low T are characteristic of a 4D gas in free space. Since (to the best of 
our knowledge) no 4D gas has been observed previously, the predicted low T and critical 
behaviors of this system are particularly interesting. 

A number of theoretical issues need to be discussed. First, one must ask whether the 
ideal gas assumption is valid; indeed, our group has previousl y e xplored the vapor-liquid 
transition of these gases using quasi-lD models of perfect ICs |42j. A concern is therefore 
that such condensation in a somewhat disordered system would either preempt the BEC or 
alter the nature of the BEC, as is the situation with 3D superfluid 4 He. However, it may 
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well be that disorder reduces the temperature of the hypothetical condensation transition 
below that of the BEC transition, so that two separate transitions occur. 

A related question is the role of screening by the medium of the interparticle interaction. 
On the bare surface of graphite, experimental and theoretical evidence has indicated that the 
well-depth is reduced by some 10 to 20% 43]. Intuition and some calculations |41( suggest 
that the attraction would be reduced by a much greater factor for gases within a nanotubes 
array. That is consistent with the neglect of these interactions in the present paper. 

The requirements for observation of this phenomenon in the laboratory are not trivial to 
satisfy. A particularly serious concern is that equilibrium is difficult to achieve. It requires 
particles to diffuse out of one IC and into another as T changes. The slow equilibration 



is a potentially fatal problem t 
be very helpful in this regard 



rat is difficult to assess; experimental diffusion data would 
44I We note that one need not imagine that particles 
are required to move macroscopic distances; one can employ tubes of quite finite length in 
order to accelerate the process of equilibration. Evidently, one should not underestimate 
this problem, which is an important concern in much broader contexts involving nanotubes, 
such as gas storage and isotope separation. 

The problem at hand raises the venerable question of the relation between BEC and 
superfluidity (SF). J4q] We know that BEC and SF appear simultaneously at the lambda 
transition of 4 He |3l, M. However, the Kosterlitz-Thouless transition of 2D films is one 

n 

involving SF without BEC |47| . Moreover, the 3D ideal bose gas exhibits BEC without 
SF, since the Landau velocity criterion is not satisfied. The latter says that the threshold 
velocity for superflow is the minimum of the ratio of the excitation energy to the momentum 
. Our unusual nanotubes-hosted BEC state would also fail this test, implying that no 
superflow is possible. Another aspect of SF is quantized circulation, which is impossible in 
this environment. Thus, we think that SF is unlikely to accompany the transition described 
above. If, however, one were to include weak interparticle interactions, SF could occur. 

The ideal nanotube radii for the lowest energy states fall in the range 0.8 to 1 nm, which is 
an advantageous size because typical samples contain many such tubes. If, for example, only 
smaller tubes are present in a particular sample, one would have to explore the behavior 
without all of the simplifying mathematics employed here. We have yet to tackle that 
problem in general, but do note some interesting results found in one case. Suppose for 
simplicity that the transverse density of states for a given sample (consisting of N/c ICs) 
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has a form (that is quantitatively wrong but qualitatively reasonable in that it consists of a 
spread of transverse energies) given by the "uniform barrier" function: 



g{E t ) 




> 



(23) 



g(E t ) = 0, otherwise 



(24) 



The result in this case (from Eq. |7|) is 



N(E) 



2N IC L m 



1/2 {{E - £<) 1/2 0[£ - E<] -{E- £>) 1/2 6[£ - E>}) (25) 



(E> - E<) nh v 2 



This expression yields the ID limiting case of Eq. 4 when the barrier width (E > — E<) falls 
to zero. In the case of nonzero width, this expression yields 3D behavior for N(E) at low 
E and ID behavior at high E. The reason for the 3D limit is that the constant transverse 
density of states in the finite uniform barrier model mimics the constant density of states of 
a 2D gas. In this model, one derives a BEC transition from this form of N(E), with thermal 
behavior quite different from that described elsewhere in this paper, i.e., when the lowest 
possible transverse energies (near E m i n ) play the dominant role in determining the behavior. 
The transition temperature of a system of N p molecules, obtained for the uniform barrier 
model with the conventional 3D theory, satisfies 



4D proportionality to iV*' 2 in Eq. El 

One might ask whether similar BEC phenomena occur in a related environment, the 
groove between two nanotubes (e.g. at the external surface of a bundle). Adsorption in this 
region has been much studied in theory and experiments j^ . We have therefore investigated 
the consequence of heterogeneity for that problem, with interesting results. One finding is 
that the symmetric situation (two identical radii) provides a local maximum (in 2D R space) 
of the transverse energy. Hence, the methodology and results of the present paper are not 
applicable, for the most part. The one related result is that there occurs in the groove case 
an anomalous signature of inhomogeneous broadening in the excitation spectrum. This is 



{k B T c f = -( 
m 

to N 2 / 3 is chai 
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qualitatively analogous to the behavior of N(E) here, except that spectroscopy is the tool 



of choice for its investigation. 



In closing, we emphasize an intriguing fact: this transition is a consequence of the nonuni- 
formity of the ICs within experimental samples. A perfectly uniform set of ICs would exhibit 
ID physical properties instead of the transition described here. Such a remarkable conse- 
quence of heterogeneity has precedents in other physical systems, of which the spin glass 
phenomenon is a well-known example |sol . I51I lo^ ] . 

We are grateful to Carlo Carraro, Vincent Crespi, Allan Griffin, Susana Hernandez, Paul 
Lammert, Ari Mizel, Aldo Migone, Paul Sokol and Flavio Toigo for helpful discussions. 
This research has been supported by the National Science Foundation. F.A. acknowledges 
funding from MIUR-COFIN 2001. 

APPENDIX: 4D IDEAL BOSE GAS 

We consider an ideal bose gas in four dimensions whose density of states N(E) satisfies 

N(E) = aE (E > 0) (A.l) 

For this system, we evaluate the number of particles, energy and specific heat in the usual 
way, using Eqs. 13-15. Straightforward calculations yield 

N p = a g 2 (Z) (k B T) 2 + N cond (A.2) 
where N con d is the number of bose condensed particles, and 

U = 2ag 3 (Z) (k B T) 3 . (A.3) 

Here the functions g 2 and #3 are second and third " Bose Einstein integrals" , 



9n{Z) = —i ttzt-Ft- (A-4) 



f°° x n l dx 
r(n) Jo Z~ l e x - 1 ' 

The fugacity Z(T) is determine from Eq. IA.21 Moreover, from this equation we obtain 
the critical temperature T c for a given number of particles N p , given that the maximum 
value of Z is 1 and g 2 (Z) is a monotonically increasing function, 
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From Eqs. IA.5I and IA.2I the condensate fraction results, 

/ = ^ = l-(^r) 2 . (A.6) 

Now we evaluate the specific heat. To do so we have to be aware that the fugacity depends 
on the temperature only for T > T c while it is equal to 1 for T < T c . Then we obtain 

CN Q^-4^(T>T c ) (A.7) 



N p k B g 2 (Z) g x {Z) 

= 6 ^rl(— ) 2 (T < T c ) (A.8) 

Note that is a continuous function of T, with its value at the critical temperature 
C V (T C ) = AAN p k B . 

In the these derivations we have used some properties of the functions g n {Z), 

9n(l) = C(n) (A.9) 

and 

Zg' n (Z) = g n - l (Z). (A.10) 
These, together with equation A. 2, give the following relationship, 

^ = -!»<*> (A-U) 
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a (A) 


e(K) 






H 2 1 


2.97 


42.8 


8.75 


-1159.2 


H 2 2 


3.23 


32.19 


9.95 


-1052.4 


4 He 3 


2.74 


16.2 


7.19 


-353.6 


4 He 2 


2.98 


16.9 


8.48 


-456.9 



TABLE I: Numerical values of Emir, and R m i n ^ denned in the text. The parameters of the potential 

1 Q, 2 3 and 3 3 
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FIGURE CAPTIONS 

1. Contours of the ground state transverse energy E t (R) — E min (relative to its minimum 
value) of an H 2 molecule as a function of the radii R = (R 1: R 2 , R3) of the tubes bordering 
the IC. R is measured relative to the minimum R m j ra ; the figure depicts behavior along the 
line (0, 1, —1). Beginning with the contour closest to the minimum (the closed curve), the 
contours correspond to relative energy 0.005, 0.01, 0.05, 0.1, 0.5, 1 and 5 K. 

2. Upper panel: transverse density of states g(E) for 4 He, computed from the distribution 
of ICs as described in the text. Lower panel: total density of states N(E). In both cases, 
the energy is measured relative to the minimum E min = — 456.88.K" 

3. Transverse and total densities of states, as in Fig.2, for H 2 . E min = — 1052. 97/T 

4. Specific heat of He (full curve) and H 2 (dots) as a function of T computed using the 
classical Boltzmann approach. 

5. Fugacity Z = exp(/3/i) as a function of T/T c for two different densities: 1 x 10~ 5 A~ 3 
(full curve) and 0.5 x 10~ 5 A~ 3 (dots) within an inhomogeneous nanotube array while the 
dashed curve corresponds to an ideal 4D gas having the same low energy density of states. 

6. Condensate fraction as a function of the temperature computed for H 2 in the nanotube 
system (full curve) and the 4D H 2 gas (dots) for iV=lxlO~ 5 A~ 3 . 

7. Heat capacity for H 2 as a function of the temperature in nanotubes (full curve) and 
for the 4D H 2 gas (dashed) (panel a) and for 4 He in nanotubes (full curve) and for the 4D 
4 He gas (dashed) (panel b). 

8. Density of H 2 molecules (full curve) and 4 He atoms (dotted curve) as a function of 
transition temperature T c . The short-dash (H 2 ) and long dash (He) curves indicate the 4D 
ideal gas limiting behavior. 

9. The occupied state energy distribution n(E,T) (relative to the total number N p ) 
defined in Eq. 15 is shown for H 2 at relative temperatures T/T c = 0.5 (full curve), 1 (dots), 
1.1 (short dashes) and 1.5( long dashes). The integral of the curve at 0.5 is the fraction (about 
0.7) of excited particles while the other curves integrate to 1, since then the condensate 
fraction vanishes. 
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